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^\ • Fermilab experiment E791 studied weak decays of , D'^ , and mesons produced in 

0^ ' collisions of 500 GeV/c negative pious with Pt and C targets. The experiment collected 

0^ . over 200 000 fully reconstructed charm decays. Four recent results are discussed here: (a) 

measurement of the form factor ratios V/Ai, A2/A1, and A^/A^ in D"*" K t^v^ and 
bJQ, Dt — > (j>(.'^Vi decays; (h) measurement of the difference in decay widths AF between the two 

^ ' D°/d'' mass eigenstates; (c) search for rare and forbidden D decays to dilepton final states; 

, and (d) search for a "Pentaquark," a bound state of csuud. 

00 

1 Introduction 

' Fermilab E791 is a charm hadroproduction experiment studying the weak decays of charmed 

1/^ ■ mesons and baryons. The experiment took data from September, 1991 to January, 1992, record- 

, ing over 2 x lO^'' interactions and reconstructing over 200 000 charm decays. This large sample 

I has led to numerous published results.El Four recent results are discussed here: (a) measure- 

G^ ' ment of the form factor ratios V/Ai, A2/A1, and A^/Ai in — > K*^l'^V(^ and Df (j)£~^ 



decays;il (^6) measurement of the difference in decay widths AF between the two D^/d'^ mass 
eigenstates;0 (c) search for rare and forbidden D decays;Q and (d) search for a "Pentaquark," a 
, bound state of csuud.!'! 

Qj I The E791 collaboration comprises approximately 70 physicists from 17 institutions. The 

' experiment produced charmed mesons and baryons using a tt~ beam of momentum 500 GeV/c 

incident on five thin target foils (one platinum, four carbon). The foils were separated along the 
\ beamline by approximately 1.5 cm such that most charm decays occurred in air rather than in 

^ • solid material. Immediately downstream of the target was a silicon strip vertex detector consist- 

ing of 17 planes of silicon strips oriented along the x, y, u and v directions, where u and v point 
± 20° from the vertical direction (y). Following the vertex detector was a spectrometer consist- 
ing of two large-aperture dipole magnets providing pj- kicks of 212 MeV/c and 320 MeV/c, and 
37 planes of wire drift chambers and proportional chambers. Downstream of the second magnet 
were two threshold Cerenkov counters used to discriminate among pions, kaons, and protons. 
Following the Cerenkov counters was a Pb/liquid-scintillator calorimeter used to measure the 
energy of electrons and photons, and an Fe/plastic-scintillator calorimeter used to measure the 
energy of hadrons. Downstream of the calorimeters was approximately 1.0 m of iron to range 
out any remaining hadrons, and following the iron were two stations of plastic scintillator - one 
X- measuring and one y- measuring- to identify muons. The experiment used a loose transverse 
energy trigger {Ej, ^ 3 GeV) that was almost fully efficient for charm decays. 

After events were reconstructed, those with evidence of a decay vertex separated from the 
interaction vertex were retained for further analysis. The analyses presented here selected events 
using numerous kinematic and quality criteria; the most important of these are listed in Table |l|. 

" C.B.P.F. (Brazil), Tel Aviv (Israel), CINVESTAV (Mexico), Puebla (Mexico), U.C.Santa Cruz, University 
of Cincinnati, Fermilab, Illinois Institute of Technology, Kansas State, University of Massachusetts, University of 
Mississippi, Princeton, University of South Carolina, Stanford, Tufts, University of Wisconsin, and Yale. 
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Table 1: Kinematic and quality criteria used to select events. 



Cut 


Typical value 


= l^^soc ~ ^prim ) 1 y "^sec + <^prim 

(transverse to D line-of- flight) 

(target) = {z^^^ — Z^^^^^ edge) /''"sec 

DIP = D impact parameter w/i/t primary vertex 
xLck/(d-o.f.) 

Lifetime = (z.^^ - Zp^j^^) • m^/p 


12 (1)0, Df), 20 (D+) 
< 250 MeV/c 
5.0 < 
< 40 ^m 
< 5.0 

< 3.0 ps (DO, L>+), < 5.0 ps {D+) 



An especially effective criterion for enhancing signal over background was "SDZ," the longitudi- 
nal distance between the production and decay vertices divided by the total measurement error 
in this quantity. 



2 and Form Factors 

Semileptonic decays such as — > K*^i^v^ and (pi^u^^ proceed via spectator diagrams. 

As such, all hadronic effects are parametrized by four Lorentz-invariant form factors: Ai{q'^), 
A2{q'^), A-i{q^), and V{q'^). Unfortunately, the limited size of current data samples precludes 
measurement of the dependence, and we assume this dependence to be given by a nearest pole 
dominance model: F{q^) = F(0)/(1 — '7^/"^poie)' where m^^^e ~ "^-^ GeV/c^ for the vector form 
factor V and 2.5 GeV/c^ for the axial- vector form factors A. Because Ai{q'^) appears in every 
term in the differential decay rate, we factor out ^i(O) and measure the ratios ry = y(0)/^i(0), 
r2 = ^2(0)7^1(0), and = A3(0)/Ai{0). These ratios are insensitive to the total decay rate 
and to the weak mixing matrix element V^g. 

To select — > K*^i^i'^ and — > cpi'^i'^ decays, we identify 3-track vertices in which one 
track is identified as a kaon and one track as a lepton. We cut on the "transverse" mass rrirp, 
where m'^ = {Ej^tj^+E^+prp)^ — ( Vk* l<ii^Pt ^ Pt)^- ■'■^ ^'^^^ expression, is the momentum 
of the neutrino transverse to the direction of the D as inferred by momentum balance. The mj^ 
distribution of semileptonic decays forms a Jacobian peak with an endpoint at m^), and thus we 
require that mj. lie in the range 1.6-2.0 GeV/c^ (1.7-2.1 GeV/c^) for the {Df) sample. We 
also require that either rrij^^ ~ m^* or ra^j^ ~ m^. The resulting samples contain very little 
background, and we do a maximum likelihood fit for the form factors using a likelihood function 
based on three angles: 9y, the polar angle in the K* {(p) rest frame between the ■k'^ (-^^) and 
the (Df); 0^, the polar angle in the rest frame between the and the {D'^); and x, 
the azimuthal angle in the (Df) rest frame between the K*^ {(p) and decay planes. The 
results of the fit to the K data for the e + /i samples combined are: ry = 1.87 it 0.08 it 0.07 
and r2 = 0.73 ± 0.06 ± 0.08. We measure rg = 0.04 ± 0.33 ± 0.29 from the n sample alone. 
The results for (pl'^u^ are: ry = 2.27 ± 0.35 ± 0.22 and r^ = 1.57 ± 0.25 ± 0.19. These 

results are compared with theoretical predictions in Table |2|; the errors in the measurements are 
smaller than the spread in theoretical predictions. 



3 D^-D^ Mixing and AF 

E791 has published a limit on the D^-D^ mixing rate using semileptonic K~£~^i'^ decaysEl 

and hadronic —>■ K~7r~^ and — > K~7r~^7r~-n'~^ decays-S The flavor of the or Z)'' when 
produced is determined by combining the D^/D^ with a low momentum pion to reconstruct 
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Table 2: Theoretical predictions for and r2- Some values are extrapolations from — 



to = 0. 



Group 



V 



E791 D+ 

ISGW20 

wsgI 

Ksl 

AW/GS0 
StechO 
BKS@ 
LMMS0 
ELC0 
APElil 
UKOCT^ 
BBD§ 
LANlE 



1.87 ± 0.11 



0.73 ± 0.10 



2.0 
1.4 
1.0 
2.0 
1.55 

1.99 ± 0.22 ± 0.33 
1.6 ± 0.2 
1.3 ± 0.2 
1.6 ± 0.3 

1 4+0-5 

2.2 ± 0.2 
1.78 ± 0.07 



1.3 
1.3 
1.0 
0.8 
1.06 

0.7 ± 0.16 ± 0.17 

0.4 ± 0.4 

0.6 ± 0.3 

0.7 ± 0.4 

0.9 ± 0.2 

1.2 ± 0.2 

0.68 ± 0.11 



E791 ^+ 

ISGWsE 

BKsEi 

LMMsli 



2.27 ± 0.41 



1.57 ± 0.31 



2.1 

2.00 ± 0.19 
1.65 ± 0.21 



+0.20 

-0.25 



1.3 

0.78 ± 0.08 
0.33 ± 0.33 



+0.17 
-0.13 



a D*^ D^7r~^ or D* D it decay. The semileptonic decays yield a 90% C.L. limit 
r^i^ < 0.50% [where r^;^ = T{D^ ^D^-> J)/T{D^ /)], while the hadronic decays yield a 
90% C.L. limit r^^^ < 0.85%. The latter limit assmnes no CP violation in the mixing and no CP 
violation in a doubly-Cabibbo-suppressed (DCS) amplitude which also contributes to the rate. 
However, CP violation is allowed in the interference between the mixing and DCS amplitudes. 
Since the DCS amplitude is in fact substantially larger than that expected from mixing, the 
presence of "wrong-sign" decays in the hadronic data - while a signature for mixing- is more 
easily interpreted as evidence for DCS decays. If we assume no mixing, then the numbers of 
wrong-sign decays observed in our data, corrected for acceptance, imply ratios of DCS decays 
to Cabibbo-favored decays of rg'^g = (0.68 tgj^ ± 0.07)% and r^^^'^ = (0.25 tg-J^ ± 0.03)%. 

Since r^-^-^ = (1/2) [(Am/F)^ + (Ar/2r)^], where Am and AT are the differences between 
the masses and decay widths of the D^/D^ mass eigenstates, the upper limit for r^^^j^ implies an 
upper limit for the difference in widths: |Ar| < 0.48 ps~^. E791 has made a direct measurement 
of AF using — > K~K~^ and — > K~-k^ decays. Since the former results in a CP- 
even eigenstate, only the CP-even component contributes and the lifetime distribution is 
proportional to e~^i*. The K^tt^ final state, however, is a CP admixture and the lifetime 
distribution is proportional to exp [- (El +r2)t/2] cosh(Art/2).y Over the range of lifetimes 
for which the experiment has sensitivity, cosh(Ari/2) ~ 1 and thus: Ej^^ — E^^ = E^ — {T^ + 
E2)/2 = (E, - E2)/2 = AE/2. 

Our samples of K^K^ and K^tt^ are shown in Eig. |l|a. We bin these events by 

reduced proper lifetime, which is the distance traveled by the candidate beyond that required 
to survive our selection criteria, multiplied by mass and divided by momentum. Eor each bin 
of reduced lifetime we fit the mass distribution for the number of signal events. Plotting this 
number (corrected for acceptance) as a function of reduced lifetime gives the distributions shown 
in Eig. Pd. Fitting these distributions to exponential functions yields F^j^ = 2.441 it 0.068 ps^^, 
F^^ = 2.420 ± 0.019 ps-\ and thus AF = 0.04 ± 0.14 ± 0.05 ps'^ This implies -0.20 < 
AF < 0.28 ps~^ at 90% C.L., which is more stringent than the constraint resulting from r^-^^. 
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Fiffure 1: D° -> K-7V+ and D° K+ mass distributions (left), and reduced proper lifetime distributions 

(right). The right-most peak in the lower left plot is due to misidentified D*^ — > K~tt^ decays. 



4 Rare and Forbidden D Decays 

E791 has searched for rare and forbidden dilepton decays of the , Df, and D^. The decay 
modes can be classified as follows: 

1. flavor-changing neutral current decays £'^i~ and ^(^s) ~^ h'^£'^i~, in which /i is a 
pion or kaon; 

2. lepton-flavor violating decays I?" — > ji^e^, D'^d^) h'^^'^e^, and D^^s) ~^ h~^~^e^, in 
which the leptons belong to different generations; and 

3. lepton-number violating decays -D^ h~i~^£~^, in which the leptons belong to the same 
generation but have the same sign charge. 

Decay modes belonging to (1) occur within the Standard Model via higher-order diagrams, but 
the branching fractions are estimatedil to be only 10~* to 10~^. This is below the sensitivity of 
current experiments. Decay modes belonging to (2) and (3) do not conserve lepton number and 
thus are forbidden within the Standard Model. However,.^ a number of theoretical extensions 
to the Standard Model predict lepton number violationJH^ and the observation of a signal in 
these modes would indicate new physics. We have searched for 24 different rare and forbidden 
decay modes and have found no evidence for them. We therefore present upper limits on their 
branching fractions. Eight of these modes have no previously reported limits, and fourteen are 
reported with substantial improvements over previously published results. 

For this study we used a "blind" analysis technique. Before our selection criteria were 
flnalized, all events having masses within a window AM5 around the mass of the D~^, Df, or 

were masked so that the presence or absence of potential signal candidates would not bias 
our choice of selection criteria. All criteria were then chosen by studying signal events generated 
by Monte Carlo simulation and background events obtained from the data. The background 
events were chosen from mass windows AMb above and below the signal window AM5. The 
criteria were chosen to maximize the ratio Ng/ y/N^, where Ng and are the numbers of 
signal and background events, respectively. Only after this procedure were events within the 
signal window unmasked. The signal windows AMs used for decay modes containing electrons 
are asymmetric around to allow for the bremsstrahlung low-energy tail. 
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We normalize the sensitivity of our search to topologically similar Cabibbo- favored decays. 
For the modes we use K^ir^ir^; for the D'^ modes we use cpn^; and for the 

we use — > K~7:~^. The upper limit on the branching fraction for decay mode X is: 

R ^norm r> ("-i \ 

— 'T7 ~ ■ "norm ) 

norm ^X 

where N^^ is the upper limit on the mean number of signal events, A^^orm number of 

normalization events, and and ^norm overall detection efficiencies. The geometric accep- 
tances and reconstruction efficiencies are found from Monte Carlo simulation, and the particle 
identification efficiencies are measured from data. 

The background consists of random combinations of tracks and vertices, and reflections from 
more copious hadronic D decays. The former is essentially flat in the reconstructed invariant 
mass, and we estimate this background by scaling the level from mass regions above and below 
the signal region AM5. The hadronic decay background in which a K is misidentifled as a 
lepton is explicitly removed via a Ktttt or KKtt invariant mass cut. The hadronic background 
in which a vr is misidentifled as a lepton cannot be removed in this manner, as the reflected 
mass and true mass are too close and such a cut would remove a substantial fraction of signal 
events. We thus estimate this background by multiplying the number of — > 7r~7r+7r"'", 
Df K~iT^iT^ , or — > vr+TT^ decays falling within the signal region AM5 by the rate for 
double particle misidentiflcation vrvr — > /i/i, /ie, or ee. The misidentiflcation rates were measured 
from data using — > K~t^~^tt~^ decays misidentifled as K~l'^l^ . Because the latter samples 
have substantial feedthrough background from the former (which is Cabibbo- favored), we do 
not attempt to establish a limit for — > K^i'^H.'^ decays. Rather, we use the observed signals 
to measure the lepton misidentiflcation rates under the assumption that all K~l'^t'^ decays 
observed arise from misidentifled i^^vr+vr"'". Most of our flnal event samples are shown in Fig. |2|, 
and all results are tabulated in Table |^. 



5 Search for the Pentaquark P^g^^^ 

E791 has searched for a "Pentaquark" P'^, which is a bound state of fii^e quarks having flavor 
quantum numbers csuud. This state was originally proposed by LipkinEd and Gignoux et alx^ 
over ten years ago, but no experimental searches have been undertaken. The is predicted to 
have a mass below the threshold for strong decay {mj^^+nip = 2.907 GeV/c^) by 10-150 MeV/c^. 
The lifetime is expected to be similar to that of the shortest-lived charm meson, 0.4-0.5 ps. We 
have searched for P^ decays into (pp-K^ and K*^K~p flnal states. 

The analysis proceeds by flrst selecting four-track vertices in which one track is identifled as 
a proton and two opposite-sign tracks are identifled as kaons. We require that either ~ 
or nij^^ ~ ''^K* remove events in which either m^^ ~ or the (j) or p momentum projects 
back to the production vertex. We normalize the sensitivity of the search to — > cj) vr^ and 
K* decays; these are topologically similar to P^ — > (j)pTT~ and P^ K*^K~p 
(except for the proton) and several systematic errors cancel. After all selection criteria are 
applied, we observe no excess of events above background in either decay channel. We thus 
obtain upper limits for the product of production cross section and branching fraction, relative 
to that for the D'^ . The expression used is (here for (f)p'ir~): 



(2) 



where Np^^^^ is the upper limit on the mean number of P^ (j)piT decays, and Njj is the 
number of events observed in the normalization channel. All numbers and the resulting limits 
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m(KV,V) 



n(K*e*e") 




m(KV.V) 



n(/x^^ ) m(e^e ) 

hvariant Mass (GeV/c^) 



n{/i*e") 



Figure 2; Final event samples for (rows 1-2), Df (rows 3-4), and D'^ (row 5) decays. The dotted curves show 
signal shape for a number of events equal to the 90% C.L. limit. The solid curves show estimated background. 



are listed in Table ^ When calculating acceptance, we assume the lifetime to be 0.4 ps. 
The limits are given for two possible values of mpo; the difference in the limits is due mainly 
to the difference in the numbers of events observed in the mass spectrum around these mass 
values. Our upper limits are 2-4% of that for the corresponding Df decay, which is similar to 
the theoretical estimate (~ 1%). 



6 Summary 

We have presented four recent results from Fermilab experiment E791: a measurement of the 
form factors governing — > K*^i^u^ and (pl'^u^ decays; a measurement of the difference 

in decay widths AT between the two mass eigenstates of D^/D^; new limits on two dozen rare 
and forbidden dilepton decays of D^, D^, and Df; and a limit on a • B for a "Pentaquark" 
relative to that for the Df. Almost all measurements and limits are superior to previously 
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Table 3: 90% C.L. upper limits for dilepton decays of D'^ /D'^ /D'^ . The right-most colmTin lists previous results. 



Mode 


Background 


Observed 


UL (AT) 


UL(S) 
(xlO^) 


1 1 

PDGQSEfl 
( X ID j 


r-\-t- 

— 




2.67 


2 


3.35 


< 


1.5 


< 1.8 


D+ - 




0.90 


1 


3.53 


< 


5.2 


< 6.6 


D+ - 




0.78 


1 


3.64 


< 


3.4 


< 12 


— 


— -I- -i- 

> IT l-L fJ- 


1.53 


1 


2.92 


< 


1.7 


< 8.7 


D+ - 




0.45 


2 


5.60 


< 


9.6 


< 11 


D+ - 




0.39 


1 


4.05 


< 


5.0 


< 11 


— 




2.40 


3 


5.07 


< 


4.4 


< 9.7 


D+ - 


> K+e+e- 


0.09 


4 


8.72 


< 


20 


< 20 


D+ - 




0.08 


1 


4.34 


< 


6.8 


< 13 


r-i + 

DJ - 


+ — 


2.00 





1.32 


< 


14 


< 59 


Dt- 




0.85 


2 


5.77 


< 


160 


— 


Dt- 




1.10 


1 


3.57 


< 


63 


- 


DJ - 


> K fi^fi^ 


1.04 





1.68 


< 


18 


< 59 


Dt- 


> i^-e+e+ 


0.39 





2.22 


< 


63 


— 


Dt- 


► K-fi+e+ 


1.15 


1 


3.53 


< 


68 


- 


Dt- 




1.65 


1 


3.02 


< 


14 


< 43 


Dt- 




0.83 





1.85 


< 


27 


— 


Dt- 




0.72 


2 


6.01 


< 


61 


— 


Dt- 


► 7r^/i+/_i+ 


1.16 





1.60 


< 


8.2 


< 43 


Dt- 




0.42 


1 


4.44 


< 


69 




Dt- 




0.36 


3 


8.21 


< 


73 








2.46 


2 


3.51 


< 


0.52 


< 0.41 




e+e— 


2.04 





1.26 


< 


0.62 


< 1.3 






2.88 


2 


3.09 


< 


0.81 


< 1.9 



Table 4: 90% C.L. upper limits on the product of cross section and branching fraction {a ■ B) for P° — > ^pvr 
decays (left table) and K*^K~p decays (right table), relative to that for —^ ^vr^ and 





P" mass (GeV/c^) 
2.79 2.87 




3.1 7.5 
0.47 0.64 
293 ± 18 


Upper limit: 


0.022 0.040 





P" mass (GeV/c^) 
2.77 2.85 


^p/^D^- 


6.1 3.5 
0.23 0.31 

725 ± 88 


Upper limit: 


0.036 0.016 
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published results. In the case of the and eight of the rare and forbidden dilepton decays, our 
limits are the first such limits reported. 
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